Sudden jumps and plateaus in the quench dynamics of a Bloch state 
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Take a one-dimensional tight binding chain with the periodic boundary condition and put a 
particle in some Bloch state, then quench it by suddenly changing the potential of some site. In 
the ensuing time evolution, the probability density of the wave function at an arbitrary site jumps 
constantly between plateaus. This phenomenon adds to another one that the survival probability 
of the particle in the initial Bloch state shows cusps periodically, which was found previously in 
the same scenario (Zhang and Yang, arXiv:1601.03569). The two, one in the real space and one 
in the momentum space, complement each other to provide a comprehensive picture of the quench 
dynamics of a Bloch state. Underlying the cusps and jumps is the exactly solvable, nonanalytic 
dynamics of a fictitious model, based on which, the locations of the jumps and the heights of the 


plateaus are accurately predicted. 


PACS numbers: 03.65.-W,02.30.Rz 


I. INTRODUCTION 


Singularities can show up in the time evolution of 
a physical quantity. A classic example is provided by 
a model solved rigorously by Stey and Gibberd four 
decades ago [1]. The model is in the Lee-Friedrichs class 
[2, 3], and is actually a paradigm for quantum decay [4— 
6]. It consists of an equaldistant quasi-continuum ex- 
tending from —oo to +00, and an extra discrete level, 
which couples to all the continuum levels with the same 
strength. Initialize the system on the discrete level and 
let it decay into the quasi-continuum. It turns out that 
the survival probability of the initial state shows cusp- 
s periodically in time. More recently, in the surge of 
nonequilibrium dynamics of many-body systems [7-11], 
Heyl et al. have discovered non-analyticities at critical 
times in the time evolution of some Loschmidt echo re- 
lated quantities [12]. In contrast to the former example, 
the non-analyticities were sought deliberately by noting 
the formal similarity between the Loschmidt amplitude 
(U;|e~*#"|U,) and the partition function Tr(e~°”), and 
by borrowing notions from equilibrium phase transitions. 
The singularities were then legitimately called dynamical 
quantum phase transitions. Thereafter, many systems 
demonstrating this novel type of phase transition in the 
time domain were found [13-18]. 

Personally, we have also found some nonsmooth dy- 
namics recently in a very simple scenario [19]. The set- 
ting is a one-dimensional tight binding chain with the 
periodic boundary condition, which is arguably the sim- 
plest model in solid state physics [20]. Initially a particle 
is put in some Bloch state with an arbitrary momentum. 
Then suddenly the potential of some site is changed. In 
the subsequent evolution, both the survival probability 
and the reflection probability, which correspond to the 
particle remaining in its initial state and being momen- 
tum reversed, respectively, show cusps periodically. The 


cusps were explained by exactly solving an idealized mod- 
el. The mechanism is completely different from those 
functioning in the examples above. 


The singularities invite experimental verifications. Un- 
fortunately, so far this goal is yet to be fulfilled. For the 
model of Stey and Gibberd, a clean realization is to put 
a two-level atom in a multi-mode cavity and to position 
it at an appropriate point [21-24]. However, to the best 
of our knowledge, this has never been attempted experi- 
mentally, possibly due to the challenge of the necessity of 
precisely locating the atom. An alternative scheme is to 
use the one-dimensional tight binding model [5]. But the 
drawback is that it is confined to the perturbative regime. 
As for the dynamical quantum phase transition, the es- 
sential difficulty is that the central quantity, namely the 
Loschmidt amplitude for a many-body system, cannot be 
measured directly. Therefore, the phase transition has to 
be inferred from some conventional quantities which are 
experimentally accessible. Some proposals in this spirit 
do exist [13, 14, 18] but experimental realization is still 
missing. 

Motivated by this problem, we have reinvestigated the 
singular dynamics of the Bloch state [19]. From the ex- 
perimental point of view, to verify the singular dynamics 
by measuring directly the two probabilities (or popula- 
tions) is not necessarily an easy task. In many circum- 
stances, it is much easier to measure local quantities. We 
have thus investigated the time evolution of the probabil- 
ity density of the wave function at an arbitrary site. From 
the theoretical point of view, this study is also worthy as 
it provides a real space picture of the dynamics, com- 
plementary to the momentum space picture established 
previously. It turns out that the finding is as interesting 
as in the momentum space. On a large time scale, the 
trajectory of the local density is characterized by plateaus 
and sudden switches between them. The lengths of the 
plateaus depend on the site. In particular, the heights 


of the plateaus are very sensitive to the location of the 
site; they vary significantly even between adjacent sites. 
Therefore, we believe measuring the local density is a 
much more feasible means to verify the singular dynam- 
ics of the particle. 

The rest of the paper is organized as follows. First, 
the cusps and plateaus will be presented side by side in 
Sec. II, demonstrating that they are two aspects of the 
same scenario. Then in Sec. III, the plateaus will be ex- 
plained by invoking the formalism developed previously 
in Ref. [19] for accounting for the cusps. After that, in 
Sec. IV, we discuss an experimental system promising for 
realizing the scenario. Finally, we conclude and discuss 
the implications of the findings in Sec. V. 


II. SUDDEN JUMPS AND PLATEAUS 


The setting is a one-dimensional tight binding chain 
with the periodic boundary condition. The original 
Hamiltonian is (A = 1 throughout this paper, and the 
hopping strength between two adjacent sites is set as the 
unit of energy) 


N- 


Ho =- X (0+1 + (+ 14), (1) 


l=0 


=. 


where |/) denotes the Wannier function at site l. By 
the periodic boundary condition, |l + N) = |l). The so- 
called Bloch states are simultaneously eigenstates of Ho 
and eigenstates of the translation operator T|) = |!+ 1). 
They are indexed by the integer k, and have the explicit 
expression of 
(k) = et! (2) 
VN , 
in the real space. Here q = 27k/N is the so-called wave 
vector. Apparently, |k} = |k + N}. The corresponding 
eigenenergies are €(q) = —2 cos q. Like the Wannier func- 
tions |l) are the basis vectors in the real space, the Bloch 
states |k} are the basis vectors in the momentum space. 
Now the scenario is as follows. Initially a particle is 
put in some arbitrary Bloch state |k;), i.e., |Y (0)) = |k}. 
Then at t = 0, the potential of some site is changed 
to U suddenly and held fixed afterwards. Without loss 
of generality, the site is assumed to be the l = 0 one. 
The final Hamiltonian which will govern the evolution 
of the wave function |W(t)) of the particle is then Hy = 
Ay + fh, with H, = U|0)(0|. The characteristic feature 
of the newly introduced perturbation is that it couples 
two arbitrary Bloch states with the same strength, i.e., 


g =U/N = (ki | Hike) (3) 
regardless of the values of kı,2, as can be easily checked 
by using the explicit expression of (2). 


The question is how does the wave function |U(t)) e- 
volve. An intuitive picture is that the particle will be 


reflected back and forth between the two groups of Bloch 


states centered at | + ;) and | — ki). Therefore, in a 
previous work [19], the two quantities 
Pit) = (+e WO), 4a) 
P(t) = (ks [Y(A))P, (4b) 


which are called the survival probability and the reflec- 
tion probability, respectively, were studied. The finding 
is that, both P; and P, show cusps periodically. Ex- 
amples are shown in Figs. l(a) and 1(c). In Ref. [19], 
the cusps and even the whole trajectories of P;,, were ex- 
plained by identifying and solving an ideal model behind. 

Quantities like P;, are defined in terms of the momen- 
tum space. They are motivated by the Rabi oscillation 
picture above. From another point of view, the newly in- 
troduced barrier at site l = 0 will excite scattering waves, 
which will propagate both forwards and backwards away 
from the barrier. Because the group velocity of a wave 
packet on the lattice chain is upper bounded, it will take 
finite times for the scattering waves to reach a site. Once 
the scattering waves come back to the barrier, they will 
generate secondary scattering waves, and so on. There- 
fore, the wave function |(¢t)) might exhibit rich temporal 
and spatial structures. 

We have thus investigated numerically the time evolu- 
tion of the probability density 


Dy(t) = (UE = [Wilt |? (5) 


of the wave function at an arbitrary site l. Correspond- 
ing to Figs. 1(a) and 1(c), we have Figs. 1(b) and 1(d), 
respectively, where the plateaus and the fast switch be- 
tween them are the most prominent feature. They are 
not rigorously plateaus—there are fluctuations. But in 
each interval, the curve fluctuates apparently around a 
horizontal line. In Figs. 1(a) and 1(b), where] = 5 < 
N = 401, it seems that the sudden jumps happen simul- 
taneously with the cusps. However, this is not the case in 
Figs. l(c) and 1(d), where 1/N œ 1/4. Further numerical 
examples indicate that the times when the jumps hap- 
pen depend on the location of the site under observation. 
This is in line with the scattering wave picture above. 


II. EXPLANATION 


Behind the plateaus and sudden jumps is actually the 
same ideal model which is responsible for the cusps [19]. 
This ideal model was motivated and justified by an im- 
portant numerical observation. That is, for a wide range 
of parameters, only those Bloch states with energies in 
the vicinity of the energy of the initial Bloch state |k;) 
participate significantly in the dynamics. This means t- 
wo groups of Bloch states, with wave vectors either in the 
vicinity of +q; or —q. Here q; = 27k; /N is the wave vec- 
tor of the initial Bloch state. Figure 2 demonstrates this 
point very well. There, for a lattice of size N = 401, and 
an initial Bloch state with k; = 100, the population on 
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FIG. 1. (Color online) Time evolution of the survival probability P; (solid lines), the reflection probability P, (dashed lines), 
and the local probability density |W;|? [see Eqs. (4) and (5) for definition]. The left two panels correspond to one case, while the 
right two panels to another case. The common parameters are N = 401, k; = 100, while the rest parameters are (U,1) = (2.5, 5) 


in left two panels, and (U,1) = (2,100) in the right two panels. 
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FIG. 2. (Color online) Time evolution of the population on 
those Bloch states outside the two slots of |k — ki| < W and 
|k + ki| < W [see Eq. (6) for definition]. The parameters are 
the same as in the left two panels of Fig. 1, ie., N = 401, 
ki = 100, and U = 2.5. 


those Bloch states outside of the two slots of |k—k;| < W 
and |k + ki| < W is studied. That is, the quantity 


do kwe? (6) 


|k—kj|>W 
|k+k;|>W 


py? (t) = 


out 


is followed. Here W is the cutoff. It is seen that even for 
W = 4, the two narrow slots, which cover less than 3% 
of the Brillouin zone, almost exhaust the probability. 
Therefore, it is tempting to construct a fictitious mod- 
el, which resembles the realistic model H f in the two 
slots but can deviate from it (even significantly) outside. 
Around the two points of +q;, the dispersion relation of 
the unperturbed Hamiltonian (1) can be linearized. The 
wave vectors are equally spaced with a gap of 6 = 2r /N 
and the eigenenergies are equally spaced with a gap of 
A = «'(q)d = 4rsinqi/N. Hence, we truncate the re- 
alistic spectrum of Ho by retaining only the two slots 
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FIG. 3. (Color online) The red solid curve depicts the disper- 
sion relation e(q) = —2 cos q of the unperturbed Hamiltonian 
(1). The red solid dots on the curve indicate the Bloch states 
(2). The blue hollow dots, which form two linear branches, 
indicate the levels in the ideal model. Note that while for the 
realistic model (1), the number of levels is finite and the value 
of the wave vector q is confined to the interval [—7z, +r], in 
the ideal model, the number of levels is infinite and q extends 
from —oo to +00. 


and then extend them into two branches of linear spec- 
tra {|R,)} and {|Ln)}, -co < n < +00. See Fig. 3. For 
n in the vicinity of zero, these states correspond to the 
Bloch states as 


|n) © lki +n), |En) © |- ki — n). (7) 
The basis states {| Rn} } and {|Ln)} are eigenstates of the 
fictitious unperturbed Hamiltonian Ho, 

Ho|Rn) =nA|Rn), HolLn)=nAlLn). (8) 
Resembling (3), the fictitious perturbation Hı couples 
two arbitrary states with the same strength 

g= (Rn, |Hi|Rns) = (En, |Hi|Lna) 
= (Rn, [Hi|Lna)- (9) 


FIG. 4. (Color online) Graphs of the functions (a) 8(z), which is piecewise linear and (b) yi(s), which are piecewise constant. 


See Eqs. (16)-(19) for definition. 


Equations (8) and (9) define the fictitious model com- 
pletely. The idea of its construction is that, the exact 
details of the Hamiltonian beyond the two slots do not 
matter much, as the population on those levels is negli- 
gible. Therefore, it is free to rectify that part to make 
analytic calculation possible. This will be justified if the 
prediction of this fictitious model agrees with the exact 
results. 

It should be stressed that although in the original prob- 
lem, there are only a finite number of (N actually) levels, 
here in the fictitious model, we have two branches of lev- 
els, each consisting of infinite levels. Moreover, in the 
original model, Bloch states |k) and |k + N) are identi- 
fied as the same state, while here 


|Rni) A |En) Vm, ne, (10a) 
|Rny) # | Fires) if nı £ N2, (10b) 
|[Lm) Z |En) if ni £ no. (10c) 


Because eventually we are to calculate the amplitude of 
the wave function at an arbitrary site l, we have to postu- 
late the expressions of the basis states {|R,,)} and {|Ln}} 
in the real space. This is done simply by generalizing (7) 
to all n, and taking (2), i.e., 


VN 


In Ref. [19], the dynamics of the initial state |Y(0)) = 
n E 


(IRn) = — exp(i2r (ki + n)l/N) = (I|Ln)*. (11) 


|k;) = | Ro) under the control of the fictitious Hamiltoni- 


an Ho + Hy has been solved. In terms of the states with 
definite parities, |A*) = (|Rn) + |Ln))/V2, the state at 
time t is in the form of 


a E 
= aldo) + Yo Yn(t)|At). 


n=—Cco 


|Y(t)) (12) 


At this point we have to introduce a very important quan- 
tity, i.e., the so-called Heisenberg time T = 27/A. This 
is the time which sets the time scale of the dynamics in 
question. It is actually the period between the cusps in 
Figs. 1(a) and l(c). For rT < t < (r+ 1)T, where r isa 
nonnegative integer, the value of wo is (s = t — rT) 


i2gs Liro 
t)=(1- ir 1 

woi)=(1- ew, 0 

and the value of Yn (n Æ 0) is 

2g e7inAs =A] aoe 

n = Fale DR p e . 14 
Palt) (1 +igT)A ( n ) e oe 
Here e~” = (1—igT)/(1+igT) is a phase factor. In (13), 


we see that Wo is a continuous but nonsmooth function 
of time t. Its trajectory on the complex plane is like 
the trajectory of a ball inside a circular billiard. This 
explains the cusps in Figs. 1(a) and l(c), as P;, = |1 + 
wo|?/4. Substituting (11), (13), and (14) into (12), we 
have 


, œO b inAs _ 1 , o i N 
NUY =a il —ird jl g e iqil+inôl —igqil—inol 
VN(I|U(t)) = ising! + e cos + (+ igh) A 2 z (e +e ) 
Ee ea eT ee ee a OA =À 
isingl +e cos qil + Gtigh (« (B(6 s) — B(6l) s) 


= ising +e”? cos ql + ame (eitla (8) + e~t4ly_(s)) | . 


In the second line we have defined the function 


ije 25 nme 


n=1 


(16) 


+e~*4!(B(—61 — As) — B(—61) as))| 


(15) 


m 


It is apparently periodic with a period of 27. In Ref. [25], 
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FIG. 5. (Color online) Comparison between the analytic predictions (red dashed lines) by (15) and the numerical exact results 
(blue solid lines) for the local probability density. The common parameters are N = 401 and U = 1. The rest parameters 
are shown in each panel. Panels (a) and (b) demonstrate that even between two adjacent sites, the evolution trajectories of 
the local probability densities can be very different. Panel (d) shows that even for large times, when the plateaus are not 
well shaped, the analytic predictions still capture the trend of evolution of \w,|? very well. Also note how the location of the 
observed site, namely the value of l, influences the lengths of the plateaus by comparing the upper panels with the lower ones. 


it is calculated to be [see Fig. 4(a) for its graph] 


B(z) = ig — 2r|z/2r| =r), z A 2ar, (17) 


0, z=2nrr, 


where |-| is the floor function. In the third line, we have 
defined the functions 


y4(s) = (+61 — As) — (+5!) — As, 
+-(s) = B(-él — As) — B(—61) — As, 


which are defined on the interval [0,7]. Using (17), it is 
easy to determine [see Fig. 4(b) for their graphs] 


(18a) 
(18b) 


0 O0<s<s 
=4 ° 19 
1) Se <s<T, ae) 
0 0<s<T-s 
aC a = 19b 
=e) oa T-s.<s<T, (13b) 


where se = 61/A. The physical meaning of se is the time 
needed for a wave packet, whose wave vector is centered 
at qi, to go from the barrier to the site l in the forward 
direction. The reason is simply that the group velocity 
of the wave packet is exactly v = e'(q;) = A/8. Simi- 
larly, the Heisenberg time T = 27/A = N/v is just the 
time needed for the wave packet to complete the loop 
and return to the barrier. The time T — se in (19b) is 


then simply the time needed for a wave packet with wave 
vector —q; to go from the barrier to the site l. 


Equation (15) is our central result. We see that the 
amplitude Ų;(t) jumps at rT + se. Between the jumps, 
the amplitude is a constant. Now we can compare the 
predictions of (15) and the numerical exact results. This 
is done in Fig. 5. We see that the analytic predictions 
based on the fictitious model (red dashed lines) agree 
with the numerical results (blue solid lines) very well. 
Generally, not only the times when the jumps occur, but 
also the heights of the plateaus, are accurately predicted 
by (15). An important prediction of (15) is that, the 
value of |Ņ;(t)|? is very sensitive to the site index l, and 
this is vividly demonstrated in Figs. 5(a) and 5(b). 


We should note the regularity of the time evolution of 
|W,|?. The expression (15) contains two parts. The first 
part, ising;l, is time independent. The second part is 
time dependent, but it factorizes into a part containing 
solely r and a part containing solely s. Hence, U;(rT +s) 
as a function of r is of the form a + be~*”®, where a and 
b are constants. Its modulus squared is then a sinusoidal 
function of r. This harmonic oscillation feature is easily 
perceived in Figs. 1 and 5. In particular, when 6/7 is a 
rational number, as in Fig. 1(d), where 0 = 7/2, |W)(rT+ 
s)|? is a periodic function of r. 


The times when the jumps happen are in accord with 
the scattering picture below. As illustrated by Fig. 6, 
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FIG. 6. (Color online) Nine equal-distant snapshots of the 
probability distribution of the wave function in the time in- 
terval [0, T]. For clarity, from bottom to top, the jth curve is 
up shifted by 2(j—1). The parameters are N = 401, k; = 100, 
and U = 1. The wave vector of the initial Bloch state is then 
qi = 2rki/N ~ w/2. The dashed line indicates the motion 
of the back scattering wave front with a constant velocity of 
v = 2sing; ~ 2. The forward scattering wave front is also 
clearly visible. 


the barrier generates forward scattering waves with wave 
vectors q ~ +q,, and back scattering waves with wave 
vectors q ~ —q;. Their wave fronts move at a constant 
velocity of 


Oe A : 
ere = 5 = 2sing (20) 


forwards or backwards. It takes them exactly the time 


N N 2r T 91 

v Aj A gu 
to return to the barrier. Then the scattering waves will 
generate secondary scattering waves which also travel at 
the velocities +v along the lattice. The times when the 
scattering waves, or the secondary ones, pass by the site 
l, are just rT se, the times when the jumps happen. 

In Fig. 6, snapshots of the probability distribution of 

the wave function in the first period [0, T] are shown. In 
the four curves at the bottom, the steep forward moving 
wave front is clearly visible. For a fixed site l, before 
the wave front reaches it, the value of |W;|? is almost 
stationary at 1/N; when the wave front crosses it, the 
value of |W;|? experiences a swift change. 


IV. EXPERIMENTAL REALIZATION 


Like the van Hove singularity in the density of states 
of a solid [26-28], here the singularity in the dynamics is 
also worth experimental verifications. 


A promising approach to observing the sudden jumps 
and plateaus is to use coupled optical waveguides [29]. 
Because of the quantum-optical correspondence, many 
coherent quantum phenomena have been successfully 
simulated using photonic guiding structures carefully en- 
gineered. In particular, some interesting phenomena with 
the tight binding model as the setting, such as the dy- 
namical localization [30, 31] and the Bloch oscillation 
[32-34], have been realized. 

In our case, the periodic boundary condition can be 
implemented by arranging the N waveguides in a circu- 
lar loop, and the initial Bloch state can be realized by 
carefully engineering the phases of the input laser beam- 
s, while the defect potential can be achieved by changing 
the refractive index of one of the waveguides, which can 
be done in several different ways [29]. The primary po- 
tential difficulty might come from the large time scale 
involved, which is linearly proportional to the lattice size 
N. In contrast, the time scales of the dynamical local- 
ization and the Bloch oscillation are independent of N. 


V. CONCLUSIONS AND DISCUSSIONS 


We have reinvestigated the single-site quench dynam- 
ics of a Bloch state in a one-dimensional tight binding 
lattice, which was previously studied by us in Ref. [19]. 
The motivation was on the one hand to get a complemen- 
tary picture of the dynamics from another perspective, 
and on the other hand to find a quantity which is more 
convenient for measurement. 

We have thus tracked the evolution of the probability 
density of the wave function at an arbitrary site, which is 
probably the quantity most readily measurable. It turn- 
s out to jump constantly, from plateau to plateau. In 
other words, its time development is not monotonic, but 
well structured. The times when the jumps happen, the 
durations as well as the heights of the plateaus, can all 
be accurately predicted by a fictitious model. Due to the 
nonlinearity of the realistic spectrum, and the finiteness 
of the number of levels in the realistic model, there are 
details beyond the predictions of the fictitious model, but 
the overall trend is well guided by its analytic prediction- 
s. The dependence of the times of the jumps on the site 
under consideration supports the scattering wave picture. 
Namely, the probability density jumps each time when a 
scattering wave passes by. 

From the point of view of thermalization or equilibra- 
tion, the model in question never thermalizes or equili- 
brates. First, in the figures, we can see that each time 
after the jump, the amplitude of fluctuation shrinks with 
time—it is equilibrating. However, this process is only to 
be interrupted by the next jump. The processes of jump 
and equilibration then alternate. Second, it is observed 
that even between two neighboring sites, the probability 
densities have completely different trajectories and never 
merge together. This means that the wave function will 
never become “even” across the lattice. Hence, the real 


space evidences corroborate the opinion that the dynam- 
ics of the model is very regular and non-chaotic. 
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